A Feshbach resonance in d-wave collisions 
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We analyse a narrow Feshbach resonance with ultra-cold chromium atoms colliding in d-wave. 
The resonance is made possible by dipole-dipole interactions, which couple an incoming I — 2 col- 
lision channel with a bound molecular state with I = 0. We find that three-body losses associated 
to this resonance increase with temperature, and that the loss feature width as a function of mag- 
netic field also increases linearly with temperature. The analysis of our experimental data shows 
that the Feshbach coupling is small compared both to the temperature and to the density limited 
lifetime of the resonant bound molecular state. One consequence is that the three body losse rate is 
proportionnal to the square of the number of atoms, and that we can directly relate the amplitude 
of the losses to the Feshbach coupling parameter. We compare our measurement to a calculation 
of the coupling between the collisionnal channel and the molecular bound state by dipole-dipole 
interactions, and find a good agreement, with no adjustable parameter. An analysis of the loss 
lineshape is also performed, which enables to precisely measure the position of the resonance. 

PACS numbers: 34.50.-s, 03.65.Nk, 67.85.-d 



I INTRODUCTION 

The study of atomic coUisional properties in ultracold 
gases has been made possible by the recent development 
of laser cooling and trapping techniques. In a Feshbach 
resonance the scattering cross section is resonantly en- 
hanced when the energy of a bound state in a given 
molecular potential matches the molecular dissociation 
limit of another molecular potential, provided both po- 
tentials are coupled. For polarized bosons at low tem- 
peratures, collision processes are dominated by s-wave 
scattering, and a Feshbach resonance results in a modi- 
fication of the s-wave scattering length. This feature is 
particularly relevant for Bose-Einstein condensates, since 
it allows for the production of BEG with tunable interac- 
tions, with for examplepossible application to interfero- 
metric measurements Feshbach resonances can also 
be used to produce molecules, by ramping the magnetic 
field through the resonance In the case of mixtures 
of fermions, s-wave Feshbach resonances are extensively 
used to study the BEC-BCS crossover around the unitary 
limit p-wave Feshbach resonances were also studied 
in the case of polarized fermions 0]. In that case, both 
the incoming open channel and the molecular bound state 
have a p-wave character. A d-wave shape resonance was 
also studied by colliding two Rb condensates at large ve- 
locities 0]. Resonances with incoming I larger than zero 
may offer a means to study the effect of anisotropic in- 
teractions in quantum degenerate gases. 

Here, we study the resonant coupling of an incoming 
channel ^ = 2 to a bound state of partial wave I — 0. 
The coupling is due to dipole-dipole interactions. As the 
bound state is localized at relatively short internuclear 
distance (typically 100 ao), the particles in the incoming 
channel need to tunnel through the centrifugal barrier 



associated to Z = 2. We shall see that tunnelling through 
this barrier has profound consequences on the tempera- 
ture dependence of the resonance width, on the temper- 
ature dependence of the three-body loss parameter, and 
even on the dynamics of atom losses. 

The paper is organized as follows. We first report on 
an experimental study of a Feshbach resonance in d-wave 
collisions between ultra cold ^^Cr atoms. We will re- 
fer to such resonance as a 'd-wave Feshbach resonance'. 
Losses are analysed, and we show that three-body pro- 
cesses lead to losses proportional to the square of the 
number of atoms, N'^, not to as is usually the case. 
We study the temperature dependence of both the am- 
plitude of the loss feature and of its width as a function 
of the magnetic field B. Both vary linearly with T. We 
then develop a theoretical model taking into account that 
the Feshbach coupling is small compared to the inverse 
(density limited) lifetime of the bound state. We finally 
compare the loss amplitude with a calculation taking into 
account the coupling of the incoming d-wave channel to 
the s-wave molecular state by dipole-dipole interactions. 
We find a very good agreement between our data and 
this calculation involving no adjustable parameter. 



II EXPERIMENTAL RESULTS 

One manifestation of the strength of dipole-dipole in- 
teractions in chromium is the existence of d-wave Fesh- 
bach resonances. In the one we study, first reported in [§], 
the incoming channel jS* = 6, Ms — — 2, mi = 1) 
is coupled to \S — 6, Ms = —5, 1 = 0,mi — 0) through 
dipole-dipole interactions. The d-wave scattering 
cross section, usually vanishingly small at low tem- 
perature because the atoms need to tunnel through 
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FIG. 1: Left: number of atoms remaining in the trap after a hold 
time of 10 s, in the vicinity of the Feshbach resonance, at temperature 
15 /iK. Right: typical time evolution of the number of atoms after the 
magnetic field is set at the magnetic field corresponding to maximum 
losses. The sold line is the result of a fit using eq. l3l . 



a centrifugal barrier, is resonantly enhanced when 
the energy of the last vibrational bound state of 
15 = 6, Ms = —5,1 = 0,mi — Q) matches the molecular 
dissociation limit of jS* = 6, Ms — —6, 1 — 2, m; = 1). 

We proceed by measuring the atom loss rate in an opti- 
cally trapped Chromium atom cloud as a function of the 
magnetic field B. We first load from a MOT up to 5 mil- 
lion ground state ^^Cr atoms at a temperature of 100 fiK 
in an optical dipole trap formed by a retroreflected far 
red detuned 35 W laser beam, following an experimental 
procedure described in [S] . To prevent inelastic collisions 
due to dipolar relaxation, we polarize the atoms in the 
lowest energy Zeeman state ms = —3, which is done by 
optically pumping the atoms with a circularly polarized 
laser at 427 nm {i.e. resonant with the ''Sa — >^ P3 tran- 
sition) in a 2.3 G homogeneous magnetic field. The next 
step is to form a dimple by adding a strong longitudinal 
confinement to the trap. For this we transfer power from 
the horizontal beam to a crossed vertical beam by turn- 
ing a half wave plate on a computer controlled rotation 
stage, in front of a polarizing beam splitter cube (as in 
0). During the linear rotation of this plate the depth of 
the horizontal beam trap decreases, forcing evaporation 
and thus decrasing the number of atoms and the tem- 
perature of the cloud. Determined by the rotation angle 
of the half-wave plate, the atom number is between 10^ 
and 3 x 10** and the temperature between 16 fiK and 2 
fjK. At this point we tune the static homogeneous mag- 
netic field to look at field dependent losses due to the 
d-wave Feshbach resonance around B — 8.2 G. We use 
RF spectroscopy to calibrate the magnetic field with an 
uncertainty of 2 niG, and resonant absorption imaging to 
extract the temperature and the atom number. 

Fig [H (a) shows the number of atoms remaining in the 
trap after 10 s as a function of B, at T = 15 /iK. The 
lineshape of the loss feature is clearly asymmetric, which 
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FIG. 2: ? ( see definition in text) measured as a function of the 
logarithm of the number of atoms. To better emphasize the power law 
dependence in Log(N)^ the different values of ^ obtained at different 
temperatures have been shifted so that they coincide for a given atom 
number. Then all values align on a linear curve of slope p—1 — 1 (solid 
line; the dotted line is a linear curve of slope p — 3). 



is in contrast to what is observed in most s-wave Feshbach 
resonances analyzed to date, but analogous to what is 
reported in the case of p-wave Feshbach resonances [4j]. 
This paper provides a full experimental and theoretical 
analysis of these losses. 

As a first step to characterize the losses, we measure 
the time evolution of the number of atoms after the mag- 
netic field is set at the value corresponding to the peak 
losses. As shown in Fig[T}3, we observe a non-exponential 
decay of the number of atoms as a function of time. As 
the temperature of the cloud does not change much dur- 
ing the decay (less than 15 percent), at any given mo- 
ment, the density is approximately proportional to the 
number of atoms; a non-exponential atom number decay 
is then clearly related to a density-dependent loss mech- 
anism. We describe the density dependent losses by the 
following rate equation: 



dt ^ 



(1) 



where n is the atom density. Fi is the one-body loss rate 
due to backgound gas collisions, and Kp represents the 
density dependent loss rate parameter. Although atoms 
are in the lowest state of energy, so that two-body losses 
are energetically forbidden, and three-body losses are the 
most likely, we do not assume that p = i. When inte- 
grated over space assuming a thermal gaussian distribu- 
tion, Eq. ([T]) leads to 
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Fi = -PpNP-^ (2) 
where /3p = p3i2(^^y2^l^3/2y-i . with (a;o,yo,zo) the 1/e 
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FIG. 3: Loss parameter K2 as a function of the magnetic field, at 
five different temperatures. Solid lines result from fits using cq. Jol . 



FIG. 4: Full Width Half Maximum (FWHM) of the loss parameter 
feature (sec Fig|3j, as a function of temperature. The solid line results 
from a fit of slope 1 to the experimental data. 



radii of the trapped gas. It is therefore interesting to plot 
^ = Log[jj^ + Ti) as a function of Log{N), as the 
slope gives a measure of p — 1. Such a plot is represented 
in FigO Fi = l/20s^^ is deduced from the time constant 
of decay when the number of atoms is very small. For 
this figure, we used data obtained at five dilTerent tem- 
peratures. As we are only interested in the power number 
dependence {p— 1) of the loss process, we have shifted the 
experimental curves vertically by a value £,o{T), so that 
all values of ^ — S,o(T) coincide for a given atom number. 
It can be seen in Fig [2] that all the data points lie close 
to the same straight line, corresponding, independent of 
the temperature, to a slope (p — 1) = 1. 

The striking conclusion of this first study is that, al- 
though two-body losses are completely excluded based on 
energy conservation arguments [3], the losses we observe, 
most likely due to three-body recombination, are well de- 
scribed by a two-body loss parameter. Such a feature will 
be explained in the next part of the paper. 

To go further in the analysis, we performed systematic 
measurements of the "two-body loss parameter" K2 de- 
scribing the three-body losses (thereafter called simply 
"loss parameter" in order to avoid confusion). 

We use two different methods to measure K2. To 
rapidly determine K2 as a function of temperature and 
magnetic field, we do not systematically record the atom 
number decay as a function of time. Rather, we use the 
number of atoms remaining after a given hold time t. 
Indeed, considering a fixed one-body loss parameter Fi, 
the number of atoms after a time t univoquely determines 
K2, as: 
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;i - exp{-Tit)) 



(3) 



The main advantage of this method is to speed up the 
data aquisition process, as a single data point provides a 
measurement of K2, instead of many data points corre- 
sponding to the time-evolution of the number of atoms. 
Results are represented in Fig [31 where (as in the rest of 
this paper) magnetic fields B are expressed in units of fre- 
quency, hv{B) = gjfisB^ where h is Planck's constant, 
gj = 2 is the Lande factor of chromium atoms in the ^83 
state, and the Bohr magneton. As evidenced in Fig[3l 
we observe a reduction of both the width and the ampli- 
tude of the loss feature when the temperature decreases, 
in spite of the increase of the phase space density. Below 
T = 2 iiK we were not able to observe any resonant loss 
feature around 8.2 G. The value of the magnetic field at 
which losses are maximal shifts with temperature, but 
the value of the magnetic field above which losses start 
to occur does not. This value is referred to as i?res, and, 
as shown below, it is the value of the magnetic field at 
which the energy of the bound molecular level exactly 
coincides with the molecular dissociation limit, i.e. the 
exact position of the Feshbach resonance. 

We also represent in Fig [1] the variation of the width 
of the observed loss feature as a function of the tempera- 
ture. We find that this width increases linearly with the 
temperature. 

For a precise determination of the magnitude of K2, 
we use a different approach: we systematically record 
the number of atoms as a function of time, and fit the 
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FIG. 5: Experimentally determined two-body loss parameter. The 
solid line results from a fit of slope 1 to the experimental data. Inset: 
Corresponding loss parameter assuming a three-body loss parameter 
n = —Ks X n^. 



data with eq. ([3]) . We report in Fig [5] the temperature 
variation of the loss parameter if measured, for each 
temperature, at the magnetic field at which losses are 
maximal. We find that within experimental uncertainty 
j^max jg proportional to T. This contrasts with the situ- 
ation of broad s-wave Feshbach resonances, for which the 
loss parameter does not depend on temperature. 

Let us now summarize the main experimental features 
of the losses close to the d-wave Feshbach resonance: (i) 
three-body losses are described by a two-body loss pa- 
rameter K2 (Fig [2]); (ii) the lineshape of the loss feature 
is asymmetric (Fig [3]); (iii) the width of the loss features 
increases linearly with T (Fig |4]) ; (iv) if™"^ varies lin- 
early with T (Fig [5]). To account for these features, we 
have developped the theoretical model described in the 
next part of this paper. 



following reaction: 



Cr + Cr 



Cr* 



(4) 



This Feshbach resonance increases the probability of 
presence of two atoms at short interatomic distance. A 
third atomic body collides with the pair of atoms be- 
fore it tunnels back out, which triggers exoergic collisions 
producing one more deeply bound molecule Cr'^'^^^'- 



Cr 



Cr 



^deep 



Cr 



(5) 



In the end, three chromium atoms (one dimer and one 
atom) have disappeared from the trap. In this frame- 
work, three-body losses can be described by the following 
Breit-Wigner cross-section: 



a{k) 



fc2 (e - £0)2/^2 + (r„,(e) + rd(n))V4 

(6) 

Here, k is the collision wavevector, e the collision en- 
ergy, and Td(n) = Pdn is the inverse life time associ- 
ated with the inelastic collision with the third body, i.e. 
the relaxation width of the bound state Crj: /3d is the 
loss parameter associated with reaction ([5]). (e — eo) is 
the energy separation to the Feshbach resonance: eo — 
AMsgjtiB{B—Bres), where AMs = 1 is the difference in 
magnetic number between the incoming channel and the 
molecular bound state, is the energy of the bound state 
relative to the dissociation limit of the incoming channel. 
Finally, UTmie) describes the strength of the resonant 
coupling associated with reaction (|4]) at a collision en- 
ergy e. It is interesting to underline the similarity of eq. 
([6|), used for the study of s-wave Feshbach resonances, to 
the standard expression for photo-association trap loss 
[l^ . which was also applied to collisions in higher par- 
tial waves. In the Wigner-threshold regime, when the 
collision energy is small, ?ir™(e) = Aie<^^W^/^ As 
a consequence, one expects that the loss parameter will 
depend on the incoming partial wave I. 

From eq. ([6]) , one can deduce the loss parameters using 



III THEORETICAL INTERPRETATION 

We interpret our experimental results using the theo- 
retical framework developed in 

0, El for s wave Fesh- 
bach resonances. Following the argumentation of these 
references, we consider that three-body losses happen in 
two steps. In a first step, two atoms collide in Z = 2. 
Tunneling to short internuclear distances is dramatically 
increased by the presence of a molecular bound state res- 
onant with the coUisional energy, and coupled to the in- 
coming pair of atoms. Such coupling is described by the 



-3 {n(T{k)vr) n — —K'i{k)n^ 



(7) 



with Vr the thermal relative velocity. The experiment 
does not directly probe Kz{k), but its averaged value 
K^{T) over the thermal distribution. If we assume a 
Maxwellian distribution at temperature T (see for exam- 
pie [il), 



nhQT 



\S{e,nf exp{^)de (8) 
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where Qt = {2T:^kBT /K^Y/"^ is the translational parti- 
tion function, with ^ = m/2 is the reduced mass. Com- 
pared to what is shown in [l^, eq. ([5]) includes a factor 
3 because we describe the loss of three particles, and 
is reduced by a factor of 2Z -I- 1 because only atoms in 
\l ~ 2,mi ^ I) contribute to losses. 

Very different results are obtained depending on the 
typical relative values of fiTm{(), Ti^d(n), and ksT. The 
goal of this paragraph is to pin-point which approxima- 
tion holds in our case; consistency will be checked below. 
hTm{i) depends on the collision energy e, but in the fol- 
lowing, we will assume that the Feshbach width can be 
characterized by the width hV^"'' = ?ir„(eS"°^) with 
= (/ -I- 1/2) X UbT. The reason for this choice of 
wiU appear later. If T™"^ >> Td{n), Td{n) can 
be neglected in the denominator of eq. (O, and the loss 
rate fi is proportional to rv^: K^[T) is then independent 
of n and atoms losses are proportional to . In order 
to reproduce our observations reported in Fig [21 we will 
therefore assume in the calculation that F™"^ << Td{n). 
Furthermore, if ?iF = ?i(r™"^ -|- Td{n)) is large compared 
to fcfiT, the width of the experimental loss feature should 
not depend on temperature (and the three-body loss pa- 
rameter should be proportional to T'). This result ap- 
plies for example to broad s wave Feshbach resonances. 
When ?IF is on the same order of magnitude as fcsT, 
the experimental width should increase with T, but not 
linearly. When fiT is very small, the lorentzian function 
in eq. ([8]) qualitatively acts as a Dirac function on the 
energy averaging function, and the experimental width 
varies linearly with T. Given the fact that the experi- 
mental width grows linearly with T, we therefore assume 
in the calculation that (?ir;^^"^ << hTd{n)) « ksT . 
Within these approximations, there is an analytical ex- 
pression for eq. ([8]): 



-Tm{£o)exp{-eo/kBT) 



(9) 



For physical interpretation, it is interesting to explic- 
itly write eq. (O, using Kz{T) as determined by eq. 
Then one simply finds for eg > 0: 



~a (nA^g) Tm{eQ)exp{-eo/kBT) x n (10) 



where A^b = 



^ , ^ is the thermal de Broglie wave- 

length, and a = 6^/2 is a numerical factor. Maximum 
losses occur for e™°^. We note that we explicitly find 
in eq. pH]) that the rate of association is, in our case, 
explicitly determined by the phase-space density, as also 
observed when the Feshbach resonance is wide [4I . 

We can now give a simple physical interpretation 
of eq. pO|) . We note \tpbound) the molecular bound 
state wavevector, Itp^) the incoming channel wavevec- 
tor, and Hdd the dipole operator coupling these two 



states. The Feshbach coupling is therefore defined as 
rm(e) = 27r |(V'f,o«nd|^fdd|^/'e)|^, which can also be writ- 
ten, in a trap, as 2nV'^p{e) = hT^le) where /ci(e) — ^ 
is the density of states in the trap {v denotes the trap 
vibrational states) and V is the matrix element for cou- 
pling between the bound state and one vibrational trap 
state [13]. If we rewrite the expression of the density of 
states as p{e) = where Av is the number of vibra- 
tional states in the trap covering an energy width set by 
Td, then N^'^'^ = Av (nA^^) is the number of atoms in 
the trap which are resonant with the bound molecular 
state of width F^. Therefore, eq. pO)) leads to an order 
of magnitude of the maximum loss rate: 



n 
— oc 



Z: 

h^Td 



X 



(11) 



We thus recover a Fermi golden rule: each trap vibra- 
tional level is an independant incoming channel weakly 
coupled to a state whose lifetime is so short that it be- 
haves like a continuum whose density of states is |^ 14 1 . 



IV QUANTITATIVE COMPARISON WITH 
EXPERIMENTAL DATA 



We now verify that eq. ([9]) quantitatively reproduces 
all the experimental findings described above. First, we 
find that K3 is proportional to 1/n: the experimental 
loss rate is proportional to n^, not to n^, although three 
particles are indeed needed for the loss to occur. In fact, 
pmaa; Td{n) correspouds to a situation where the 
collisionnally limited lifetime of the molecules due to re- 
action ([5]) is much shorter than the timescale for this 
molecule to be coupled back to the continuum through 
Feshbach coupling. For this reason, the experimental loss 
rate is only set by the slow reaction (|4]), which is a two- 
body process. This is why F^; does not appear in eq. 

Only in the limit of very small densities should we 
recover a loss rate proportional to n^. This result is in 
agreement with our experimental findings evidenced in 
Fig [2] three-body losses are in the present case best de- 
scribed by a two-body loss parameter K2 = x n. This 
is a signature that indeed F™"^ << Tdin): the precise 
study of the dynamics of the atom number decay there- 
fore provides an indirect signature of the short lifetime of 
the molecular bound state. Note that in the inset of Fig 
[21 we also report a three-body loss parameter to enable 
the comparison with three-body loss parameters in other 
systems. 

We also find that the experimental lineshapes (shown 
for example in Fig [3]) are well fitted by eq. ([9]). Tem- 
perature is used as a free parameter in these fits, and 
this yields fitted temperatures close to the ones exper- 
imentally determined by analysing expanding clouds in 
free fall, at any temperature between 3 pK and 15 fiK. 
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Note that the width of the theoretical lineshape is pro- 
portional to the temperature, consistent with our obser- 
vations reported in Fig [D The temperatures given by 
the fit are nevertheless systematically smaller than the 
measured temperatures by about 20 percent. Such dis- 
crepancy may arise from the fact that, when the loss rate 
is high (eo ~ £0"*^); there is a slight increase in tempera- 
ture, which in turn raises the loss parameter. This effect 
may decrease the half- width at half maximum of the loss 
feature that we measure. Another possibility that we 
cannot exclude is that our temperature measurements 
are slightly distorted by the presence of a small stray 
magnetic field with spatial curvature. 

We now turn to the temperature dependence of the 
loss parameter. As shown in Fig \5\ the maximum loss 
rate parameter K^'^^ increases linearly with tempera- 
ture. This is in agreement with eq. (|10p which shows 
that the maximal theoretical two-body loss parameter 
{i.e. at e™"^) varies as T'-i, therefore as T for a d-wave 
Feshbach resonance. 

Thus, eq. ^ reproduces all the experimental features 
with a very satisfying accuracy. We can therefore in turn 
use it to deepen the knowledge on this d-wave Feshbach 
resonance. For example, from the analysis of the loss 
parameter, we can measure precisely the position of the 
Feshbach resonance. This position is not the magnetic 
field at which losses are maximal, but rather, the position 
at which losses start to occur, which precisely happens 
when the bound molecular state crosses the open channel 
asymptote. Our measurements indicate: Bres = (8.155± 
0.01 5) G. Fields are precisely calibrated in the vicinity of 
the Feshbach resonance by rf atomic spectroscopy. The 
precision in our measurement is a factor of 6 better than 
what was previously reported [15 1. 

Comparing our experimental result to eq. Q, we also 

deduce that « (5.5 ± 1.2 ± 2.5) x 10"^ (|^)^^^- 

The error bars represent respectively the statistic error 
bar, and the systematic error bar, dominated by the un- 
certainty on the oscillation frequencies of the trap. For 
example, at a temperature of 8 /i K, r™''^/27r 10 
Hz. It is striking that, although the Feshbach coupling 
is so small that at any given time only a small fraction 
(« 10"'*) of the atoms can experience it, the Feshbach 
resonance can nevertheless lead to substantial losses, if 
one waits long enough. Thermalization constantly in- 
sures that there are atoms whose collision energy matches 

^max 
^0 

Having determined F™"^, we can now check for the 
consistency of the approximations which lead to eq. ©. 
Clearly TiF™"^ << ksT. The exact timescale for vi- 
brational relaxation is not known in the case of Cr, but 
it is expected to be on the same order of magnitude as 
for other bosonic species, i.e. rd{n)/2Tr « 1 kHz. The 
approximation hT^'^^ << KTdin) << ksT is therefore 
most likely valid. More precisely, HTdin) << ksT is 



obvious at the experimental densities, while F™"^ << 
Fd(n) implies molecular lifetimes of less than 10 ms. The 
good agreement between our theoretical model and our 
measurements thus indirectly provides an upper-bound 
for the coUisionnally limited molecular lifetime of the 
bound state. 

There are three main conclusions to this study. First, 
our analysis allows us to precisely measure the position 
of the resonance. Second, the width of the Feshbach res- 
onance Tm has been measured, and it is very small. This 
arises from the fact that, in a d wave Feshbach resonance, 
particles need to tunnel through a centrifugal barrier to 
reach the molecular bound state. We also show that 
Tm << rd(^), which indicates that it should be diffi- 
cult to create a large number of cold molecules by sweep- 
ing a magnetic field through the Feshbach resonance, as 
achieved in other experiments using a s wave Feshbach 
resonance A third information is that K2 oc T; going 
to lower and lower temperatures insures that losses can 
be strongly reduced. 

Eq. (fTO|) also shows that, in the regime where the col- 
lisional lifetime of the Feshbach molecules l/F^ is short, 
it is not likely to use such d-wave Feshbach resonances in 
a BEC to tailor anisotropic interactions: for the interac- 
tions to have a substantial effect on the trapped atoms, 
one needs to have F™"^ larger than the trapping frequen- 
cies, but then the lifetime of the cloud is shorter than the 
trap period when (nA^^) « 1. 



V THEORETICAL ESTIMATE OF THE 
FESHBACH COUPLING 

Losses as predicted by Eq. (fTU]) depend on a single 
parameter Vm{^o). We now turn to the comparison be- 
tween the measured value of r„i(eo) and a theoretical 
estimate obtained by calculating the coupling between 
the collisional channel and the bound molecular state by 
dipole-dipole interactions. As stated above, the param- 
eter F„(e) is defined as 27r| < ^bound\Hdd\'4^e > P, and 
depends on the matrix element of the dipole-dipole inter- 
action between the bound molecular state and the d-wave 
collisional state with energy e. We recall the expression 
of the dipole-dipole operator 



TT ^ 2 Mo ('5'i-5'2) 
Hdd = 4/^5 - 



i{Si.R){S2.R) 



■47r 



i?3 



(12) 



where Si and ^2 are the spin operators of the two in- 
teracting atoms with interatomic distance R along the 
internuclear axis R, R = RR. In tensor operator de- 
scription, one has 



dd 



(13) 



where the projection of the tensor rank is taken on the 
internuclear axis. By separating angular and radial parts, 
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as allowed by the Born-Oppenheimer approximation, the 
relevant states of an atomic pair can be written as 



11 >= 15 = 6, Ms 



-6, £^2, me = l> Fi{R) (14) 



for the collision channel, with Fi (R) the energy normal- 
ized radial wavefunction, and 

|2 >= \S = 6, Ms = -5, £ = 0, mi = 0> F2(R) (15) 

for the bound molecular state, which is the first bound 
level starting from the dissociation limit. F2{R) is the 
radial wavefunction of this molecular state. The angular 
part of the matrix element of eq. (fT5|) is easily evaluated 
using standard tensor operator technique, remembering 
that the angular momentum projections appearing in eqs 
(|14|15p are taken on a fixed axis (determined by the static 
magnetic field): 



< {Si S2) S Ms £ m,|{Si.S2}o|(5i ^2) S' M'g £' m', >= 

S{Ms + nii, M's + m'^){-l f-^'^ 



the energy of the first bound level of the S = 6 molecu- 
lar potential of Cr2, -23MHz, which is deduced from the 
measurements of Feshbach resonances of ref. [1] and cor- 
responds to a scattering length of 103 ao. The Cg value 
(733 atomic units) is also taken from |8i]. 

The calculation of the free and bound radial wavefunc- 
tions is performed by inward numerical integration of the 
radial differential equations, starting from large R values. 
For the bound state, asymptotic exponentially decreasing 
behavior is imposed and the energy is found iteratively 
to ensure that the solution vanishes at the node position. 
For the continuum, two solutions are calculated, with re- 
spective asymptotic behavior sin(fca;) and cos{kx); the 
wavefunction is the linear combination of these solutions 
vanishing at the node position. The square of the radial 
integrals 



He) 



FiiR)F2iR) 



i?3 



R^dR 



(17) 



^5(25* + 1)(2S" + l){2£ + 1){2£' 



Si{Si + l)i2S-_ 




gives the energy dependence of r„ 
obey the Wigner law 



(e), which is found to 



(18) 



in a large energy range (the relative discrepancy is less 
than 1% at 100/i K). For ksT = 8fj. K, the calculation 
gives 



The radial wavefunctions are described in a simple, 
purely asymptotic model. The method is based on the 
concept of nodal lines , |T3| ■ For energy values close to 
the dissociation limit (either above or below), the nodal 
structure of the radial wavefunctions in the inner part of 
the potential varies very slowly with energy. In a first ap- 
proximation, the radial position of the nodes of the wave- 
functions are located on straight lines as a function of col- 
lision energy, which characterize entirely the inner part 
of the potential, which can thus be ignored. The radial 
Shrodinger equations are solved in the asymptotic part 
only (i.e. at large interatomic distances), with the con- 
straint that the wavefunction vanishes on a given nodal 
line, chosen in the inner/ asymptotic frontier region. The 
nodal lines have to be adjusted to reproduce experimen- 
tal results (scattering length values, Feshbach resonances, 
near-threshold bound level positions, etc.). We use here 
a simplified version of the method, which consists of a 
R^^ potential limited at short range by an infinite repul- 
sive wall, having a tuneable position chosen to reproduce 
the scattering length of the system. The model depends 
thus on two parameters, the van der Waals Cg coefficient 
and the position of the wall, i.e. of the chosen nodal line 
[l8| . The node position is adjusted to roughly reproduce 



(7.3 ±0.3) X 10 



-5 



(19) 



which is very close to the experimental value, (5.5 ± 1.2 ± 
2.5) X 10^^ 



VI CONCLUSION 

In this paper, we have analysed a Feshbach resonance 
resonantly increasing three-body losses associated to col- 
lisions in d— wave. As the atoms have to tunnel through 
a centrifugal barrier before experiencing the resonant 
bound state, the Feshbach coupling is small compared 
to the inverse lifetime of the bound molecular state. As 
a consequence: 

• atom losses due to three-body recombination are 
described by a two-body loss parameter K2] 

• as the width of the Feshbach coupling is small com- 
pared to ksT, the observed resonance width in- 
creases linearly with T; 

• the loss parameter peak value (as a function of the 
magnetic field) increases linearly with T. 

All these features are well accounted for by a theo- 
retical model with no free parameter. We directly re- 
late the amplitude of the losses to the coupling between 
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the incoming channel and a molecular bound state, due 
to dipole-dipole interactions. In addition, our measure- 
ments provide a better determination of the position of 
this d-wavc Feshbach resonance, in addition to an upper 
bound to the lifetime of the molecular bound state. 
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